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We show that the dispersion properties imposed by an external periodic potential ensure both 
energy and quasi-momentum conservation such that correlated pairs of atoms can be generated by 
four wave mixing from a Bose-Einstein condensate moving in an optical lattice potential. In our 
numerical solution of the Gross-Pitaevskii equation, a condensate with initial quasi-momentum fco 
is transferred almost completely ( > 95%) into a pair of correlated atomic components with quasi- 
momenta k\ and k2, if the system is seeded with a smaller number of atoms with the appropriate 
quasi-momentum k\. 
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Bose-Einstein condensates in optical lattices provide 
flexible systems for studying the behavior of coherent 
matter in periodic potentials. Considerable attention is 
given to studies in regimes far from the region of validity 
of mean-field analysis and the Gross-Pitaevskii equation 
0, but also the highly non-trivial mean-field dynamics 
has been and continues to be subject to theoretical and 
experimental investigations 0, S 01 • 

The process we wish to consider is a four wave mixing 
(FWM) process, which transfers pairs of atoms coher- 
ently from an initial momentum state fco to new states 
with momenta fci and fc2. We consider a Bose-Einstein 
condensate in a quasi-lD geometry and we will consider 
only the longitudinal dynamics of the condensate. This 
geometry is relevant, e.g., for atomic wave guides and 
atom interferometers based on atom chips ^ . 

In Refs. 0, ) it was shown that nonlinear interaction 
originating from the s-wave scattering between atoms 
leads to depletion of the condensate and emission of pairs 
of atoms at other momenta when a continuous matter 
wave beam passes through a finite region with enhanced 
interactions. For a larger condensate, however, the pro- 
cess will not be effective unless it conserves both energy 
and momentum, i.e., the waves must be phase-matched 
over the extent of the sample. We shall show how the 
characteristic energy band structure in a one-dimensional 
optical lattice can be used to ensure both energy and 
quasi-momentum conservation, i.e., phase-matching of 
the FWM process. 

Our tailoring of the dispersion properties of matter 
waves by an external potential is inspired by approaches 
to non-linear optics, which employ various means to en- 
sure phase- matching, e.g., of the FWM process [sll^liol. 
A similar phase-matched FWM process has been used 
to explain giant amplification from semiconductor mi- 
crocavities, where the polariton dispersion properties can 
be controlled by the strong photon-exciton coupling . 
We also note that a recent analysis |tl2f| of the break-up 
of a bright matter wave soliton was analyzed in terms of 
dispersion and phase-matching. Phase-matched FWM 



has been realized in collisions of two condensates in two 
dimensions 13. I14IIir| . but in the present paper we show 
that the process can take place with atomic motion along 
a single direction, for example inside an atomic waveg- 
uide. 

The basic idea of our proposal is illustrated in Fig. ^ 
In a periodic potential V{z), the energy spectrum con- 
stitutes a band structure, and the figure shows the low- 
est energy band for the corresponding linear Schrodinger 
equation. When two atoms with momentum fcp collide 
and leave with momenta fci and fc2 momentum conserva- 
tion requires 



2fco = fci + ^2 modulo Q, 



(1) 



where Q is a reciprocal lattice vector. In the periodic po- 
tential the energy does not vary quadratically with the 
wave number, and as indicated by example in Fig. ^ it 
is possible to identify sets of wave numbers with conser- 
vation of the total energy 



2eo = £l + £2- 



(2) 



To investigate the effectiveness of this degenerate 
FWM process, we have performed a mean-field analy- 
sis of the dynamics of the condensate based on the one- 
dimensional Gross-Pitaevskii equation 



in—— 

dt 



2m dz'^ 

where the periodic potential V{z) is given by 
V{z) = -/3ERCOs{2kLz), 



(3) 



(4) 



where En = h?k\/2m is the recoil energy. The periodic 
potential can be generated with a standing wave of a laser 
with wavelength A = 27r/fci. The factor 



7 = gN/A^ 



(5) 



describes the strength of the nonlinearity, where N is the 
total number of atoms in the condensate, A±^ is the area 



2 




FIG. 1: Band structure for atomic motion in the periodic 
potential Eq. Quasi- momentum conservation and energy 
conservation is fulfilled in the crosses where two atoms with 
momentum ko collide and separate at momenta ki and k2 
illustrated in the figure. 
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of the transverse ground state, and g relates to the s- 
wave scattering length Og and the mass m of the atoms 
as g = Anh^as/m. In our calculation presented below we 
assume (3 — 1/2 and 7 = i0.8ER/kL corresponding to 
N ~ 100000 atoms confined to a transverse area of 
A± = 42^m^ and a grid of 512 periods of the potential. 
The effective area A± = :^\A + '2asN/L results from 
a gaussian variational ansatz [l^ ll^ , to the radial dis- 
tribution in a harmonic potential with ujj_ = 271 4,4, Hz 
and a constant longitudinal density N / L. 

An eigenstate of the condensate with quasi- 
momentum fco is found using the method of steepest de- 
scent in imaginary time, assuming a solution according 
to Bloch's theorem on a single period of the lattice poten- 
tial. Subsequently, a seed at variable ki is applied giving 
the following initial wave function 



l + ae'(''i"'=°)^ 



^o(^), (6) 



where a = 0.1 has been used in our calculations with the 
Gross-Pitaevskii equation. This wave function does not 
fulfill Bloch's theorem, and we hence restrict ourselves 
to values of fco and fci with interference patterns which 
are periodic on an extended grid of 512 periods of the 
potential. To test the importance of phase-matching in 
the FWM process we propagate the wave function '^init 
on this grid, with different values for the seeded momen- 
tum component ki. As a function of time, we can ob- 
serve the evolution of the Gross-Pitaevskii wave function 
and build-up of amplitude at different momenta. We are 
particularly interested in the quasi-momentum regions 
around fco, ki and k2 — 2fco — ki. Let ipikjt) denote the 
Fourier transform of the time-dependent Gross-Pitaevskii 
wave-function The distribution in momentum 

space folded into the single Brillouin zone from fc = to 



FIG. 2: Population of different momentum components (a): 
Pkg it), (b): Pk^ (t) and (c): Pk^ (t) as a function of time and as 
function of the seeding wave vector ki. The calculations are 
performed with a potential modulation /3 = 1/2, — 100000 
atoms, and an initial wave vector of ko = 0.672fc_L. When 
(ko, fci, fe) fulfill the phase- matching conditions in Eq. (0 and 
which is the case for fci — 0.289fci, 1.055fci, the conden- 
sate originally having wave vector ko is efficiently transferred 
into a set of atomic clouds with wave vectors ki and ^2. 



Q — 2kL is given by the following expression 

fk+Ak/2 
//c-Afc/2 



^fc+Afe/2 

Pk{t)^Y. |^(fc + 7ig,i)l'dfc, (7) 

„ Jk-Ak/2 



where we have introduced a sampling over a narrow mo- 
mentum window with Afc = ki^/32. 

Fig. [2Ia) shows the part of the condensate, Pka{t), 
remaining at the initial quasi-momentum fcg when the 
condensate is seeded with different values of fci. The 
most important features in the figure occur when the 
condensate is seeded with ki = 1.055A;l and ki — 
0.289fcL. The original condensate fraction at ko is al- 
most completely depleted, and strong growth of the 
population of the seeded momentum state is shown in 
Pfej (t) in part b, accompanied by simultaneous growth 
in the phase-matched component k2 — 2kQ — fci, shown 
as Pfe2 (t) in part c of the figure. Comparing the set 
{ko,ki,k2) = (0.672A:i,1.055A:i,0.289fcL) with the set of 
phase-matched wave vectors in Fig. ^ we find extremely 
good agreement and we conclude that the structure at 
ki = 1.055A;l is indeed due to the phase-matched FWM 
process. 

The remaining structures in Fig. |3 are less prominent 
but for instance the structure in Fig.[21[b) after 25 msec at 
k[ — 0.492/ci can be identified as a double FWM process 
with the following steps: 2ko k[ + k2 and ko + k'2 
k[ + k'^, where fcg = ^2 + (^0 ^ These steps do not 
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FIG. 3: Identification of momentum states (feojfcijfe) fulfill- 
ing the phase-matching condition. Crosses: sets of phase- 
matched values (fco, ki, ^2) identified by numerical solution of 
the Gross-Pitaevskii equation l^. Full line: phase-matching 
condition derived from band structure based on simple linear 
Schrodinger equation. Dashed line: phase-matching condi- 
tion when the effects of interactions have been included as in 
Eq. 0. Although the simple band structure calculations are 
not exact theories for the phase-matching condition there is 
still excellent agreement with the full mean-field calculation. 

conserve energy but the resulting six wave mixing process 
(S/cq 2k[+k'^) is resonant. A closer examination indeed 
confirms that the k[ component is twice as populated as 
the ^3 population in the region of the bright spot on the 
figm'e. It is again the phase-matching that is responsible 
for the significance of this amazing process. 

We have performed calculations with the Gross- 
Pitaevskii equation Q for various fco and identified the 
phase-matched sets of wave vectors {ko, ki, k2), plotted 
as crosses in Fig. |21 To quantitatively understand the 
occurring sets of wave vectors we performed a simple 
band structure calculation based on a linear Schrodinger 
equation and extracted the sets of wave vectors fulfilling 
energy and quasi-momentum conservation as illustrated 
in Fig. n The phase-matched quasi-momenta fco occur 
within the interval: fcL/2 < fcg < kzm, where fc^ = (5/2 
and kzm is a point of zero effective mass in the band 
structure = 0, < k^m < k^). Those sets of 

wave vectors are plotted as the full line in Fig. 13 This 
simple procedure produces results very much in agree- 
ment with those obtained by the full numerical solution. 
For completeness we present also a calculation including 
the effect of interactions. To calculate the band struc- 
ture for quasi-momentum states of atoms outside the fco 
condensate state, we assume that these are particle-like 
excitations, and we solve the following equation: 

( " £^ + ^('^ + 27|*o(z)p)«(z) = eu{z), (8) 

sometimes referred to as the Popov approximation to the 
coupled Bogoliubov-de Gennes equations for the prob- 
lem. The mean-field interaction term 27|^'o(z)P depends 
on fco and amounts to less than a five percent correction 



to V{z). Hence, we identify ^£'0(2;) and we solve Eq. © 
for each quasi-momentum fco and derive band structures, 
from which we extract the phase- matched pair fci , fc2 of fi- 
nal momenta. Sets of (fcoi ^i; fe) found with this method 
are shown by the dashed curve in Fig. |31 As expected 
the interactions only slightly change the phase-matching 
condition. 

Our original expectations were that the phase- 
matched, degenerate FWM could be achieved in a per- 
turbative regime with only few atoms expelled from the 
original condensate. When the process is phase-matched, 
however, the calculations show an extremely high con- 
version efficiency (up to 95%). Furthermore the pop- 
ulations Pki{t), displayed more clearly in Fig. ^a)-(c) 
for the set of phase-matched wave vectors (fco, fci, fc2) = 
(0.672fci, 1.055fcL, 0.289fcL), show clear oscillatory be- 
havior. Such Rabi-like oscillations are normally met in 
transitions between discrete states, but, as illustrated 
in Fig. ^ increasing fci (and decreasing fc2) lowers the 
energy of both states, whereby energy conservation re- 
stricts the coupling to a narrow part of the momentum 
state continuum, in which case the dynamics passes to 
Rabi oscillatory dynamics The expected frequency 
of the Rabi-oscillations is proportional to N. Fig. Efd)- 
(f ) illustrate Pkg (t) , Pki {t) , Pk^ {t) for the phase-matched 
(fco,fci,fc2) = (0.672fcL,1.047fcL,0.297fcL) from a Gross- 
Pitaevskii simulation with half the amount of atoms as 
compared with Fig.01a),(b),(c), and, indeed, we observe 
approximately half the Rabi-oscillation frequency. Note 
the slightly different values of (fci, fc2) due to the effect of 
the interactions on the phase-matching condition. 

Assuming the rather strict final state selectivity due 
to energy and momentum conservation we have per- 
formed a number state analysis in a few mode basis 
based on the three populated quasi-momentum states. 
We have thus written the wave function as 5* = 

S"feo,nfci ,"fc2 ''"'!0'"'=l ^"'=2 \^^o^ ^fel ; '^k2 X^ri Cn\N — S ~ 

2n, s + n,n > and calculated the evolution of the c„ co- 
efficients, starting with N — s atoms in fco and s in fci, 
and evolving under a Hamiltonian that removes pairs of 
fco atoms and creates fci , fc2 pairs. This calculation con- 
firms the oscillatory behavior between the mean number 
of atoms in the fco state and in the fci,fc2 states. While 
the numerical solution of the Gross-Pitaevskii equation 
shows that even a tiny seed leads to the same conversion, 
but delayed with respect to the results of a larger seed, 
our number state calculation shows that it is necessary 
to seed the condensate by a finite amount to avoid the 
decoherence due to strongly n-dependent coupling am- 
plitudes. With sufficient seed our number state analysis 
reproduces the macroscopic conversion and population 
oscillations observed in the mean field analysis. 

In conclusion the process of phase-matched FWM in a 
BEG in a periodic potential can be extremely efficient. 
When seeded, up to 95% of the atoms originally in the 
condensate with wave vector fco can be transferred into 
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will be dominant. 

K.M. Hilligs0e acknowledges NKT Academy for finan- 
cial support and Kevin Donovan for carefully reading the 
manuscript. 



0.025 0.05 0.025 0.05 0.025 0.05 

t[s] 

FIG. 4: The condensate evolution based on the Gross 
Pitaevskii equation for sets of phase-matched wave vectors: 
(a): Pfeo=o.672fcj:(t), (b): Pfei^i.ossfet (*) and (c): Pfc2=o.289fci, (i) 
for /3 = 1/2 and iV = 100000. The evolution of the phase- 
matched components of the condensate with half the num- 
ber of atoms (A'^ — 50000) is shown in (d): Pfco=o.672fcj^ (t), 
(e): Pfci=i.047fei (t) and (f): Pfe2=o.297fci (t). Half the oscilla- 
tion frequency is observed since the oscillation frequency is 
proportional to the number of atoms. 



the correlated pair of states with ki and ^2- The system 
can be used as a source of correlated atomic clouds where 
pairs can be easily separated due to the rather large mo- 
mentum difference between ki and ^2- Because of the 
coupling to a very narrow continuum of states the sys- 
tem performs Rabi-oscillations between the fco and the 
k\,k2 states, a behavior confirmed by the number state 
analysis with few modes. 

Potentially better descriptions than our simple Gross- 
Pitaevskii equation have been proposed to describe 
transversely confined elongated condensates. Our main 
mechanism relies on the non-trivial band structure, but 
not on its particular shape, and we therefore believe that 
our proposal should remain generally valid. Transverse 
confinement is, however, an important issue, and it is a 
natural extension of our theory to consider transverse 
excitations and, more generally, motion in 2D or 3D 
lattices, where energy conservation and phase-matching 
may lead to a range of interesting solutions. 

Theories 0,^3 liave proposed and experiments [l9ll20| 
have shown that condensates moving in periodic po- 
tentials become unstable for certain ranges of quasi- 
momenta. These results are linked with the energy and 
momentum conserving processes identified in this paper, 
but they also involve the detailed properties of the trans- 
verse confinement (171] . Our calculations assume a lower 
density of atoms than in the experiments reported in |0| , 
and in this regime we expect that only few atom pairs will 
be spontaneously scattered, and that our seeded process 
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